Mon. Not. R. Astron. Soc. OOO.ITHmIO Printed 20 January 2013 (MN WT^K style file v2.2) 



Mass accretion to young stars triggered by flaring activity 
in circumstellar disks 



i ■ Salvatore Orlando 1 * Fabio Reale 2 ' 1 , Giovanni Peres 2 ' 1 , Andrea Mignone^ 

O ' 

1 INAF - Osservatorio Astronomico di Palermo, Piazza del Parlamento 1, 90134, Palermo, Italy 
5_j ' 2 Dip. di Fisica, Universita degli Studi di Palermo, Piazza del Parlamento 1, 90134, Palermo, Italy 
O l' 3 Dip. di Fisica Generate, Universita degli Studi di Torino, via Pietro Giuria 1, 10125, Torino, Italy 



Accepted. Received 



ABSTRACT 

Young low-mass stars are characterized by ejection of collimated outflows and by cir- 
cumstellar disks which they interact with through accretion of mass. The accretion 
builds up the star to its final mass and is also believed to power the mass outflows, 
which may in turn remove the excess angular momentum from the star-disk system. 
However, although the process of mass accretion is a critical aspect of star formation, 
some of its mechanisms are still to be fully understood. A point not considered to 
date and relevant for the accretion process is the evidence of very energetic and fre- 
quent flaring events in these stars. Flares may easily perturb the stability of the disks, 
thus influencing the transport of mass and angular momentum. Here we report on 
three-dimensional magnetohydrodynamic modeling of the evolution of a flare with an 
idealized non-equilibrium initial condition occurring near the disk around a rotating 
magnetized star. The model takes into account the stellar magnetic field, the gravita- 
tional force, the viscosity of the disk, the magnetic-field-oriented thermal conduction 
(including the effects of heat flux saturation) , the radiative losses from optically thin 
plasma, and the coronal heating. We show that, during its first stage of evolution, the 
flare gives rise to a hot magnetic loop linking the disk to the star. The disk is strongly 
perturbed by the flare: disk material evaporates under the effect of the thermal con- 
duction and an overpressure wave propagates through the disk. When the overpressure 
reaches the opposite side of the disk, a funnel flow starts to develop there, accreting 
substantial disk material onto the young star from the side of the disk opposite to the 
flare. 

Key words: accretion, accretion disks - MHD - stars: circumstellar matter - stars: 
flare - stars: pre-main-sequence - X-rays: stars. 



1 INTRODUCTION 

Classical T Tauri Stars (CTTSs) are young low- mass stars 
activ e ly accreting mass from a surro unding disk (|Hartmannl 
1 19981 : iBouvier fe AppenzelleJ 120071 ). On the basi s of the 
large ly accepted magnetospheric accretion scenario |Koenigll 
1 19911 ) . the accretion process from the disk is regulated by the 
stellar magnetic field which is strong enough to disrupt the 
inner part of the disk at a distance of a few stellar radii (the 
truncation radius) where the magnetic pressure is approxi- 
mately equal to the total gas pressure. The field guides the 
circumstellar material along its flux tubes toward the central 
protostar, around free-fall velocity, terminating in a shock at 
the photosphere. This scenario is supported by much obser- 
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vational evidence, among others, the photospheric magnetic 
field of a few kG that has been detecte d in a number of 
CTT Ss by exploiting the Zeeman effect l|johns-Krull et all 
Il999l) . 

One of the fundamental issues in the magnetospheric 
accretion scenario is the mechanism of inward transport of 
matter and outward transport of angular momentum needed 
to explain the final mass and angular momentum of mature 
stars. The turbulence in the disk has been proposed to be the 
main responsible of the required outward angular momen- 
tum transpo rt, thus controlling the m ass accretion on the 
central star (|Shakura fc Sunva ev 1973). Nowadays, there is 
a large consensus in the literature that the turbulence is 
mainly driven by the magnet orotational instability (MRI; 
iBalbus fc Hawlevl Il99ll . 1 19981 ). MRI operates in regions 
where the ionization is sufficient to couple the gas to the 
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magnetic field and leads to the transfer of angular momen- 
tum along the ma gnetic field lines connecting gas located 
in di f ferent orbits jHawlev et ahl 1 1995: Br andenburg et aD 
ll995l ; lArmitagelll998l ; lHawlevll2000l ; lHirose et alj|200a ). De^ 
spite the significant theoretical progress achieved in this 
field, the details of this transport mechanism have yet 
to be fully understood, and it is not clear if MRI alone 
is able to account for the mass accretion and for the 



removal of excess angular momentum ( King et alj 120071 ; 
iFromang fc Papaloizoul 120071 ; iBodo et allhooa ). Given the 



complexity of this mechanism and our poor knowledge of 
its details, the efficiency of angular momentum transport 
within the disk is often modeled in the literature by in- 
cluding a viscosity in the disk modula ted via an analogue 
of th e Shakura-Sunyaev Q-parameter (|Shakura fc Sunvaevl 
1 19731 ) which can be expresse d in terms of the fluctuating 
velocity and magnetic field jRomanova et alj 1 20021 , 120031 : 
iKulkarni fc Romanovall200Sl ). 

On the other hand, observations in the X-ray band have 
shown that pre-main-sequence stars are strong sources with 
X-ray luminosity 3-4 ord ers of magnitude greater than that 
of th e present-day Sun (|Getman et al.1 120051 ; lAudard et al.l 
l2007h . The source of this X-ray radiation is a plasma at 
1 — 100 MK in the stellar outer atmospheres (coronae), 
heated by mag netic activity analogous to the solar one but 
much stronger (jFeigelson fc Montmerle] ri999). X-ray flares 
are violent manifestations of this magnetic activity and are 
triggered by an impulsive energy input from coronal mag- 
netic field. X-ray observations in the last decades have shown 
that flares in CTTSs have amplitudes lar ger than solar ana- 
logues and occur much more frequently (iFavata et al.l 2005; 
iGetman et ai1l2005l : lAudard et atboOTl : lAarnio et al.ll2010h . 
Examples of these flares are those collected by the Chandra 
satelli te in the Orion sta r-formation region (COUP enter- 
prise; iFavata et al.l [2005). The analysis of these flares re- 
vealed that they have peak temperatures often in excess of 
100 MK, are long-lasting, and are confined in very long mag- 
netic structures - extending for several stellar radii - which 
may connect the star's photosphere with the accretion disk 
(structures that could be of the same kind as those which 
channel the plasma in the magnetospheric accretion). 

At the present time, it is unclear where these flares oc- 
cur. The differential rotation of the disk together with the 
interaction of the disk with the magnetosphere may cause 
magnetic reconnection close to the disk's surface, triggering 
large-scale flares there. In this case, the flares may perturb 
the stability of the circumstellar disk causing, in particular, a 
strong local overpressure. The pressure gradient force might 
be able to push disk's matter out of the equatorial plane into 
funnel streams, thus providing a mechanism to drive mass 
accretion that differs from that, commonly invoked in the 
literature, based on the disk viscosity which determines the 
accumulation of disk matter and t he increase of gas pres - 
sure close to the truncation radius l|Romanova et alj [2002). 
Bright flares close to circumstellar disks may therefore have 
important implications for a number of issues such as the 
transfer of angular momentum and mass between the star 
and the disk, the powering of outflows, and the ionization 
of circumstell ar disks thus influe ncing also the efficiency of 
MRI (see also lAarnio et alj|20ich . 

In this paper, we investigate the effects of a flare on the 
stability of the circumstellar disk with a three-dimensional 



(3D) magnetohydrodynamic (MHD) simulation. We model 
the evolution of the star-disk plasma heated by a strong en- 
ergy release (with intensity comparable to that of flares typi- 
cally observed in young stars), close to a thick disk surround- 
ing a rotating magnetized CTTS. The model takes into ac- 
count all key physical processes, including the gravitational 
force, the viscosity of the disk, the magnetic-field-oriented 
thermal conduction, the radiative losses from optically thin 
plasma, and the coronal heating. The wealth of nonlinear 
physical processes governing the star-disk system and the 
evolution of the flare made this a challenging task. In Sect. 
[2] we describe the MHD model and the numerical setup; in 
Sect. [3] we describe the results; in Sect. [4] we discuss the 
implications of our results and draw our conclusions. 



2 PROBLEM DESCRIPTION AND 
NUMERICAL SETUP 

Our model describes a large flare in a rotating magne- 
tized star surrounded by a thick quasi-Keplerian disk. The 
flare occurs close to the inner portion of the disk, within 
the corotation radius (i.e. where a Keplerian orbit around 
the star has the same angular velocity as the star's sur- 
face), where accretion s treams are expected to originate 
|Romanova et al.1 |2002| ; iBessolaz et~alf 120081 ). The mag- 
netic field of the star is aligned dipole-like, with intensity 
B ~ 1 kG at the stellar surface according to observations 
(|johns-Krull et al.lll999h . The fluid is assumed to be fully 
ionized with a ratio of specific heats 7 = 5/3. 



2.1 MHD equations 

The system is described by the time-dependent MHD equa- 
tions in a 3D spherical coordinate system (R, 8, (f>), extended 
with gravitational force, viscosity of the disk, thermal con- 
duction (including the effects of heat flux saturation) , coro- 
nal heating (via a phenomenological term), and radiative 
losses from optically thin plasma. To our knowledge, this 
is the first numerical time-dependent global simulation of 
the star-disk system that takes into account simultaneously 
all key physical ingredients necessary to describe accurately 
the effects of a flare on the structure of the circumstellar 
disk. The time-dependent MHD equations written in non- 
dimensional conservative form are: 



dp 
dt 



+ V • (pit) = , 
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are the total pressure, and the total gas energy per unit 
mass (internal energy, e, kinetic energy, and magnetic en- 
ergy) respectively, t is the time, p = pmHrm is the mass 
density, p = 1.28 is the mean atomic m ass (assuming metal 
abun dances of 0.5 of the solar values; lAnders fc Grevessd 
1989), niH is the mass of the hydrogen atom, rh is the hy- 
drogen number density, u is the gas velocity, f is the viscous 
stress tensor, g = V<& g is the gravity acceleration vector, 
$ g = —GM*/R is the gravitational potential of a central 
star of mass M*, G is the gravitational constant, R is the 
distance from the gravitating center, T is the temperature, B 
is the magnetic field, F c is the thermal conductive flux, A(T) 
represents the optically thin radiative losses per unit emis- 
sion measure derived w ith the PINTofALE spectral code 
(|Kashvap fc Drake! l2000h and w ith the APED VI. 3 atomic 
line database jSmith et al.ll200ll ). assuming the same metal 
abunda nces as before (as ded uced from X-ray observations of 
CTTSs; iTeileschi et al. I l2007l ), and Q(R,0,cf>,t) is a function 
of space and time describing the phenomenological heating 
rate (see Sect. l2.3 (l. We use the ideal gas law, P — (7 — l)pe. 

The viscosity is assumed to be negligible in the extended 
stellar corona and effective only in the circumstellar disk. 
In order to make the transition between the disk and the 
corona, we track the original disk material by using a tracer 
that is passively advected in the same manner as the den- 
sity. We define Cdisk the mass fraction of the disk inside 
the computational cell. The disk material is initialized with 
Cdisk = 1, while Cdisk = in the extended corona. Then 
the viscosity works only in zones consisting of the original 
disk material by more than 99% or, in other words, where 
Cdisk > 0.99. The viscous stress tensor is defined as 



T = f?v 



(V«) + (Vw) T -|(V-u)J 



(5) 



where rj v = v v p is the dynamic viscosity, and v y is the 
kinematic viscosity. Several studies suggest that turbu- 
lent diffusion of magnetic field in the disk is determined 
by the same processes that determine turbulent viscos- 
ity, leading to angular momentum transport in the disk 
( Bisn ovatvi-Kogan fc Lovelacell200ll ). Thus, we assume that 
turbulent magnetic diffusivity is of the same order of mag- 
nitude as turbulent viscosity like in the Shakura-Sunyaev 
model l|Shakura fc Sunvaevl 19731 ). The kinematic viscosity 



is expressed as i/ v = acj /fix , where c s is the isothermal 
sound speed, Qk is the Keplerian angular velocity at a given 
location, and a < 1 is a dimensionless parameter regulat- 
ing the efficiency of angular momentum transport within 
the disk, which can be expressed in terms of the fluctuating 
velocity and magnetic field. The parameter a varies in the 
range . 01 — 0.6 in the Shakura-Sunyaev accretion model 
(Balbus] [2003l ). In our simulation, we assume a = 0.02. 

The thermal conduction is highly anisotropic due to the 
presence of the stellar magnetic field, the conductivity be- 
ing highly red uced in the direction transverse to the field 
l|Spitzedll962l ). The thermal flux therefore is locally split 
into two components, along and across the magnetic field 
lines, -F c = Fu i + F± j. The thermal conduction formulation 
also accounts for heat flux saturation. In fact, during early 
phases of flares, rapid transients, fast dy namics and steep 
therm al gradients are known to develop l|Reale fc Orlando! 
2008). Under these circumstances the conditions required 
for classical "Spitzer" heat conduction may break down to 



the exte nt that the plasma thermal conduction becomes flux 
limited (I Brown et al .11 19791 ). The two components of ther mal 
flux are therefore written as (|Orlando et alJl20Qgt |2010| ) 



Fn = 



F± 



[<?spi]|| feat] I 



1 1 

+ 



(6) 



[?spi]_L [<?satU 



to allow for a smooth transition between the classical and 
saturated conduction regime, where [g sp i]| and [q S pi]± repre- 
sent the classica l conductive fl ux along and across the mag- 
netic field lines (|Spitzerl 1962) 



[g S pi]|| = -«||[VT]| 
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"11 



(7) 



T 1 / 2 B 2 



[vn 



where [VT]|| and [VT]j_ are the thermal gradients along and 
across the magnetic field, and ku and n± (in units of erg s~ 
K _1 cm -1 ) are the thermal conduction coefficients along 
and across the magnetic field, respectively. The saturated 
flux along and across the magnet ic field lines, [g sa t] | and 
[<2sat]±, are (|Cowie fc McKedl 19771 ) 



feat] || = -sign ([VT] ||) 5ipp4, 
[?sat]± = -sign ([VT]j_) 5ippcs 



(8) 



where ip is a number of the order of unity; we set w = 1 ac- 
cording to the values suggested for stellar corona e (|Giulianil 
ll984lBorkowski et al. I ll989l ; iFadevev et afll2002j ). 

The calculation s were performed using PLUTO 
( Mign one et al.l |2007| ) . a modular, Godunov-type code for 
astrophysical plasmas. The code provides a multiphysics, 
multialgorithm modular environment particularly oriented 
toward the treatment of astrophysical flows in the pres- 
ence of discontinuities as in the case treated here. The code 
was designed to make efficient use of massive parallel com- 
puters using the message-passing interface (MPI) library 
for interprocessor communications. The MHD equations are 
solved using the MHD module available in PLUTO, con- 
figured to compute intercell fluxes with the Harten-Lax- 
Van Leer approximate Riemann solver, while second or- 
der in time is achieved using a Runge-Kutta scheme. A 
Van Leer limiter for the primitive variables is used during 
the heating release (at the very beginning of the simula- 
tion for t < 300 s) and a monotonized central difference 
limiter at other times. The evolution of the magnetic field 
is carried out adopting the constrained transport approach 
( Bals ara fc Spiced 19991 ) that maintains the solenoidal condi- 
tion ( V ■ B — 0) at machine accu racy. We adop t ed the "mag- 
netic field-splitting" techniqu e (|Tanakalll994l ; IPowell et al.l 
1 19991 ; IZanni fc Ferreiral 120091 ) by splitting the total mag- 
netic field into a contribution coming from the background 
stellar magnetic field and a deviation from this initial field. 
Then only the latter component is computed numerically. 
This approach is particularly useful when dealing with low- 
B plasma as it is the ca se in proximity of the stellar surface 
l|Zanni fc Ferreirall2009T ). PLUTO includes o ptically thin ra- 
diativ e losses in a fractional step formalism (|Mignone et al.l 
120071) . which preserves the 2 nd time accuracy, as the ad- 
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vection and source steps are at least of the 2 nd order ac- 
curate; the radiative losses A values are computed at the 
temperature of interest using a table lookup/interpolation 
method. The thermal conduction is treated separately from 
advection terms through operator splitting. In particular we 
adop ted the super-time-stepping technique l|Alexiades et al.l 
1 19961 ) which has been proved to be very effective to speed 
up explicit time-stepping schemes for parabolic problems. 
This approach is crucial when high values of plasma tem- 
perature are reached (as during flares), explicit scheme be- 
ing subject to a rather restrictive stability condition (i.e. 
At < (Ax) 2 /(2r/), where r\ is the maximum diffusion coeffi- 
cient), as the thermal conduction ti mescale T CO nd is shorter 
than the dynamical one Tdyn (e.g. Huieirat fc Camenzindl 
l2000l ; lHuieirat|[2fJ05l ; lOrlando et al]|2005l . l2008h ; in particu- 
lar, during the early phases of the flare evolution, we find 



Tcond / Tdyn 



10" 



The viscosity is solved with an explicit 



scheme, using a second-order finite difference approximation 
for the dissipative fluxes. 



2.2 Initial and boundary conditions 

A star of mass M» = 0.8Mq and radius R* — 2Rq is lo- 
cated at the origin of the 3D spherical coordinate system 
(R,9,<j>), with the rotation axis coincident with the normal 
to the disk midplane. The rotation period of the star is as- 
sumed to be 9.2 days. The initial unperturbed stellar atmo- 
sphere is approximately in equilibrium and consists of three 
components: the stellar magnetosphere, the extended stellar 
corona, and the quasi-Keplerian disk. 

The pre-flare magnetosphere is assumed to be force- free, 
with dipole topology and magnetic moment /xb aligned with 
the rotation axis of the star. Thus the two field components 
in spherical coordinates are 



2/xb cos 9 



fiB sm f 



(9) 



R 3 R 3 ' 

The magnetic moment is chosen in order to have a magnetic 
field strength of the orde r of 1 kG at the stella r surface 
according to observations l|johns-Krull et aT1ll999l ). 

The initial corona and disk are set in order to satisfy me- 
chanical equilibrium involving centrifugal, gravitational, and 
pressure gradient forces. In particular, we adopted the ini- 
tial conditions introduced bv lRomanova et al l l|2002f) (where 
the reader is referred to for more details) and describing the 
star-disk system in quiescent configuration. These conditions 
assume that initially the plasma is barotropic and that the 
disk and the corona are both isothermal with temperatures 
Td and T c , respectively. With these assumptions, the ther- 
mal pressure at any point of the spatial domain is given by 



P = 



Po exp 
Po exp 



2T C 



2T d 



P^ Po 
P^Pq 



(10) 



where Po = 2p c kB.T c / (pmn) is the initial pressure at the 
boundary between the disk and the corona, p c is the mass 
density of the corona close to the disk truncation radius Pd, 
fee is the Boltzmann constant, T is the function 



T = {k - 1 



GM, 

Rd 



($g+<*>c) 



(11) 



that the disk is slightly no n-Keplerian (we set k = 1.01 
as in lRomanova et aUl2002r ), $ c is the centrifugal potential 
written as 



k 



GM, 

Rd 

GM, 



1 + 



P; 



2R\ 



r^R d 
Rd 



(12) 



and r = Psin# is the cylindrical radius. The mass density 
depends on the pressure given in Eq. [10] at any point of the 
spatial domain 

[imnP 

P(P) = ' 



2T C 
(imH-P 



2T d 



P<Po 
P>Po 



(13) 



The angular velocity of the plasma is given by 



GM, 

J 'd 



1/2 



r^Rd 



f,GM,\ 1/2 ^ n 



(14) 



k is a constant of the order of 1 that takes into account 



In our simulation, we assume the isothermal disk to be 
cold (T d = 8 x 10 3 K), dense (n d ranges between 5 x 10 
and 4 x 10 12 cm -3 ), and to rotate with angular velocity 
close to the Keplerian value Ok- The rotation axis of the 
disk (coincident with the rotation axis of the star) is aligned 
with the magnetic moment [ib- The disk is initially trun- 
cated by the stellar magnetosphere at the radius P d where 
the total gas pressure of the disk equals the magnetic pres- 
sure, p + pu 2 — B 2 /8ir; for the parameters adopted here, 
Rd = 2.86 R,. In our simulation, the corotation radius 
Pco = (GM,/n 2 ,) 1/3 = 9.2 R,, where fi„ is the stellar angu- 
lar velocity. The corona is initially isothermal with T c = 4 
MK and at low densitjQ with n c ranging between « 10 8 and 
10 9 cm -3 . As shown in Eg. 1141 we allow the corona to be ini- 
tially rotating with angular velocity equal to the Keplerian 
rotation rate of the disk in order to have approximately equi- 
librium conditions and reduce the effects of transients caused 
by the i nitial differential rotat i on between the disk an d the 
corona (|Romanova et ai1l2002l ; IZanni fc Ferreirall2009h . 

Figure [1] shows the initial condition together with the 
numerical grid adopted in our simulation. The computa- 
tional domain extends between P m i n = R, (i.e. the inner 
boundary coincides with the stellar surface) and P max 
I 1 It , in the radial direction, and encompasses an angular 
sector going from Oxnin = 5° to # m ax = 175° in the angu- 
lar coordinate 9, and from m i n = 0° to ^ ma x = 360° in 
the angular coordinate <f>. The inner and outer boundaries 
in 9 do not coincide with the rotation axis of the star-disk 
system to avoid extremely small 8<f) values, vastly increasing 
the computational cost. On the other hand, all the evolution 
relevant for this study never involve portions of the domain 
close to the star-disk rotation axis. 

The radial coordinate R has been discretized on a log- 
arithmic grid with the mesh size increasing with R (see 
Fig. [l|, giving a higher spatial resolution closer to the star 
as it is appropriate for simulations of accretio n flows to a 
star with a dipole field (|Romanova et ai1l2002l ). The radial 



1 Note however that the density of the outer corona is not im- 
portant for the flare evolution. 
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grid is made of ]Vr = 80 points with a maximum resolution 
of A_R = 4.8 x 10 9 cm close to the star and a minimum 
resolution of AR = 6.4 x 10 10 cm close to the outer bound- 
ary. The angular coordinate 9 has been discretized uniformly 
with Ng = 90 points, giving a resolution of AO — 2°. The 
angular coordinate (j> is nonuniform with the highest reso- 
lution in an angular sector of 180° placed where the flaring 
loop and the stream evolve (see bottom panel in Fig. [T]). 
The </>-grid is made of jV^ = 110 points with a maximum 
resolution of A<j) = 2° and a minimum resolution of 9°. The 
numerical grid is not static but tracks the hot loop and the 
stream as the calculation progresses, in such a way that the 
loop and the stream evolve in the portion of the domain with 
the highest spatial resolution (namely that with A(j> = 2°). 

The boundary conditions at the stellar surface 
i? m in amount to assuming that the infalling material 
passes through the surface of the star as done by 
Romanov a et al.l 1 2002!') (ou t flows boundary condition; se e 
also Romanova et al. | jfloj ; iKulkarni fc Romanova! l2008h , 
thus ignoring the dynamics of the plasma after it impacts on 
the star. Zero-gradient boundary conditions are assumed at 
the outer boundary of the coordinate R (R msx ) and at the 
boundaries of the angular coordinate 9 (0 m in and # m ax). Fi- 
nally, periodic boundary conditions are assumed for angular 
coordinate 6 (</> m i n and ma x). 



2.3 Coronal heating and flare 

The phenomenological heating is prescribed as a component, 
describing the stationary coronal heating, plus a transient 
component, triggering a flare. The former component works 
at temperature T 4 MK and is chosen to balance exactly 
the local radiative losses, leading to a quasi-stationary ex- 
tended corona during the whole simulation. The transient 
component describes the injection of a heating pulse at the 
surface of the accretion disk. The flare could be driven by 
a stressed field configuration resulting from the twisting of 
magnetic field lines induced by the differential rota tion of the 
inner rim of the disk and the stellar photosphere (|Shu et al.l 
1 19971 ). The heat pulse has a 3D Gaussian spatial distribu- 
tion located at the disk border at a distance of 5 R, (namely 
well below the corotation radius R co = 9.2 _R») with width 
a — 2 x 10 10 cm. Its intensity per unit volume is Ho = 32 
ergs cm -3 s _1 . The pulse starts at the beginning of the sim- 
ulation (t — 0) and is switched off completely after 300 s. 
The flare parameters are analogous to those adopted in a 
one-dimensional (ID) hydrodynamic model that reproduces 
the evolution of a large flare observ ed in an Orion young 
star from COUP (jFavata et all 120051 ). The total energy re- 
leased during our simulated flare is En = 10 36 ergs, namely 
the same order of magnitude of the ene rgy involved in the 
brightest X-ray fl ares observed in COUP (|Favata et al ]|2005l ; 
IWolk et alj|2005h . 
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3 RESULTS 

3.1 Evolution of the flare and X-ray emission 

We followed the evolution of the star-disk system for ~ 2 
days, focusing on its effects on the disk structure. Figure [2] 
shows the distributions of density and temperature in (7?, z) 



Figure 1. Initial conditions and computational domain. Upper 
panel: volume rendering of the mass density, in log scale, at the 
beginning of the simulation (t = 0). The selected magnetic field 
lines marked in red describe the initial dipolar magnetic field. 
The yellow sphere at the center of the spatial domain represents 
the central protostar. Middle panel: slice in the (x, z) plane of 
the mass density distribution with overplotted the computational 
grid. Bottom panel: as in the middle panel for the slice in the 
equatorial plane (x,y). 
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Figure 2. Close view of the flaring loop during its early evolution. The figure shows slices in the (R, z) plane passing through the middle 
of the flaring loop, reporting the distributions of density (upper panels) and temperature (lower panels). The slices encompass an angular 
sector going from the rotation axis to the disk midplane. The white lines represent sampled poloidal magnetic field lines. The arrows 
represent the poloidal flow velocities. The magenta line delimits the region with plasma < 1 (on the left of each panel). 



slices passing through the middle of the heating pulse dur- 
ing its evolution. During this impulsive phase, the local mag- 
netic field is perturbed by the flare and an MHD shock wave 
develops in the magnetosphere above the disk and propa- 
gates radially away from the central protostar in regions 
where f3 > 1. At the same time, the sudden heat deposition 
determines a local increase of temperature (up to a maxi- 
mum of 800 MK) and pressure (above 6000 dyn cm -2 ). The 
dense disk material is heated and expands in the magneto- 
sphere with a strong evaporation front at speeds above 4000 
km s _1 . The fast thermal front propagates toward the star 
along the magnetic field lines and reaches the stellar surface 
on a timescale of « 1 hour (left panels in Fig. [2]). At this 
point, a hot magnetic tube (loop) of length L ~ 10 12 cm 
(i.e. co mparable to loop len gths inferred from COUP obser- 
vations; |Fa\2jtoirtj2[2005|) is formed, linking the inner part 
of the disk to the star's photosphere. The loop is illustrated 
in the left panels of Fig. [3] showing a cutaway view of the 
star-disk system at t — 1.2 hours (upper panel; the flaring 
loop is marked in red), and a schematic view of the system 
during the evolution of the flaring loop (lower panel) . Due to 
the efficient thermal conduction and radiation, the plasma 
begins to cool immediately after the heat pulse is over, and 
the maximum temperature rapidly decreases to ~ 50 MK in 
ss 1 hour. The disk evaporation is fed by thermal conduction 
from the outer hot plasma, even after the end of the heat 
pulse. 



Figure [2] also shows that the overheating of the disk 
surface at the loop footpoint makes a significant amount 
of material expand and be ejected in the magnetosphere. 
A small fraction of this evaporated disk material streams 
along the loop accelerated toward the star by its gravity, 
fills the whole tube in ~ 10 hours, and the density keeps 
increasing throughout the loop up to values ranging between 
10 9 and 10 10 cm 3 at its apex. On the other hand, most of 
the evaporated disk material is not efficiently confined by the 
magnetic field and channeled into the hot loop but is rather 
ejected away from the central star in the outer stellar corona, 
carrying away mass and angular momentum (see Fig. [2]). 
Due to the high values of ft there, the magnetic field lines 
are dragged away. The outflowing plasma is supersonic and 
its speed is of the order of the local Alfven speed ma ~ 300 
km s _1 . We note that our results on the dynamics of the 
outflowing plasma a re similar to those fo und with an MHD 
model proposed by Ha vashi et al.l (1 19961 ) to describe hard 
X-ray flares in protostars observed by the ASCA satellite. 

From the model results, we synthesized the X-ray emis- 
sion originating from the flare, applying a methodology anal- 
ogous to that described in the literature in the context of 
the s t udy o f novae and supernova remnants (|Orlando et al.l 
l2ood l2o'ogh . In particular, we first calculate the emission 
measure in the j-th domain cell as emj = n^Vj (where 
nfjj is the hydrogen number density in the cell, Vj is the 
cell volume, and we assume fully ionized plasma). From 
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Figure 3. Evolution of the bright flare and accretion stream. The upper panels show cutaway views of the star-disk system showing 
the mass density of the disk (green) at t = 1.2 hours (on the left) and, from the opposite side, at t = 36 hours (on the right) since the 
injection of the heat pulse. The upper left panel also overplots the three-dimensional volume rendering of the plasma temperature (in 
MK; see color table in the upper left corner of the panel), showing the flaring loop (in red) linking the inner part of the disk with the 
star. The upper right panel shows the accretion stream triggered by the flare in the side of the disk opposite to the flaring loop. Selected 
magnetic field lines are overplotted in red. Lower panels show schematic views of the system during the evolution of the flaring loop (on 
the left) and in the subsequent period, during the evolution of the stream (on the right). 



the values of emission measure and temperature in the cell, 
we synthesize the corresp onding X-ray spectrum, using the 
PINTofALE spectral code (|Kashvap fc Drakel200ol ) with the 
APED VI. 3 atomic line database, and assuming the same 
metal abundances of the simulation, namely 0.5 of the so- 
lar values, as deduce d from X-ray observations of CTTSs 
|Telleschi et al.ll2007fl . We integrate the X-ray spectra from 
the cells in the whole spatial domain and, then, derive the X- 
ray luminosity by integrating the spectrum in the [0.6 — 12] 
keV band. 

Figure U shows the evolution of the maximum temper- 
ature and emission measure of the flaring loop, and its X- 
ray lightcurve. The peak temperature of the flare is reached 
very soon at t ~ 100 s, the emission measure peaks later 
at t » 1000 s. The X-ray luminosity evolves as the emission 
measure, reaching a peak value of Lx ~ 6.5 x 10 32 ergs s _1 at 
t « 500 s, namely after the end of the heat pulse. Then the 
luminosity (and the emission measure) decreases by more 
than two orders of magnitude until t ~ 3 x 10 3 s, is steady 



till t ~ 2 x 10 4 s, and decreases again afterward. The peak 
X-ray luminosity of the simulated flare is consistent with 
the values derived for t he brightest X-ray flares observed in 
COUP (see Table 1 in lFavata et al.ll2005l ) which range be- 
tween 10 32 ergs s" 1 (source COUP 752) and 8x 10 32 ergs s" 1 
(source COUP 1568). On the other hand, it turns out that 
3/4 of the COUP flares have a peak X-ray luminosity ap- 
proximately one order of magnitude lower than that sim- 
ulated here although the total energy released during the 
simulated flare (namely En = 10 36 ergs, see Sect. 12. 3[) is 
comparable with the median total energy of the fla re in- 
ferred from the COUP observations (|Wolk et alj|2005l ). This 
may be due to the fact that the simulated flare lasts for a 
time interval significantly shorter than those typical of stel- 
lar flares and is only partially confined by the magnetic field 
(see discussion below). 

The evolution of our simulated flare has significant dif- 
ferences with re spect to the evol ution of flares simulated 
with ID models (|Reale et al.| [T988). In fact, at variance with 
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Figure 4. Upper panel: evolution of maximum temperature (red 
line) and emission measure (blue line) of the flaring loop. Lower 
panel: X-ray lightcurve of the flare in the [0.6 — 12] keV band. 



ID models where the flare is assumed to be fully confined 
by the magnetic field, in our simulation the flare is only par- 
tially confined by the magnetic field at the loop footpoint 
anchored at the disk surface. There, a substantial amount 
of the evaporated disk material escapes in the outer stel- 
lar magnetosphere and does not fill the post-flare loop. We 
conclude therefore that our results could be intermediate 
between those fou nd with models of fully confined flares 
1988) and those of models of unconfined flares 
20021 ) which show that the the flare evolution 

On 



Rcale etal 



Reale et al 



is much faster than that observed in the confined case, 
the other hand, for the purposes of this work (namely the 
study of the perturbation of the disk by a bright flare), we 
have used a simplified and idealized configuration of the ini- 
tial stellar magnetic field (namely a magnetic dipole), whilst 
many observations indicate that the stellar atmospheres are 
permeated by magnetic fields w ith a high degree of com- 
plexity (e.g. iGregorv et "all |2010| and references therein). In 
particular, the magnetic field in proximity of a heat release 
(due to magnetic reconnection) near the disk is expected to 
be more complex than that modeled here. In the presence of 
complex magnetic field configurations, the magnetic struc- 
ture hosting the flaring plasma is expected to confine more 
efficiently the hot plasma, producing a flare evolution more 
similar to that described by ID models. 



It is worth emphasizing that our 3D simulation is fo- 
cused on the effects of the flare on the stability of the disk 
and does not pretend to describe accurately the evolution of 
the flaring loop. Nevertheless, we show here that, even if the 
flare evolution is not described accurately, the length, maxi- 
mum temperature, and peak X-ray luminosity of the flaring 
loop reproduced by our simulation resemble those derived 
from the analysis of t he brightest X-ray flares observed in 
young low-mass stars l|Favata et al.l|2005l ). We are confident, 
therefore, that the flare simulated here is appropriate to in- 
vestigate the effects of bright flares observed in young stars 
on the stability of the disk. 



3.2 Dynamics of the accretion stream 

During the flare evolution, the injected heat pulse produces 
an overpressure in the disk at the footpoint of the loop. This 
overpressure travels through the disk and reaches the oppo- 
site boundary after ~ 5 hours, where it pushes the plasma 
out to form an intense funnel stream. Figure [5] shows the 
distributions of density and pressure in (R, z) slices passing 
through the middle of the stream. The overpressure wave 
triggering the stream is evident in the bottom panels. This 
new intense stream flows along the dipolar magnetic field 
lines and impacts onto the stellar surface « 25 hours after 
the injection of the heat pulse. The right panels in Fig. [3] 
show a cutaway view of the star-disk system (upper panel) 
after the impact of the stream onto the stellar surface and 
a schematic view of the system during the stream evolution 
(lower panel). As a result, the stream accretes substantial 
mass onto the young star from the side of the disk opposite 
to the post-flare loop. Our 3D simulation follows the evo- 
lution of the accretion stream for additional 23 hours, for 
a total of 48 hours. In this time lapse the stream gradually 
approaches a quasi-stationary condition. 

We analyzed the dynamics of the stream by deriving 
the forces at work along a fiducial magnetic field line nested 
within the stream. Initially the stream is triggered by a 
strong pressure gradient due to the overpressure wave orig- 
inating from the flare. The pressure gradient force drives 
the material out of the disk and channels it into a funnel 
flow. Then the gravitational force accelerates the escaped 
material toward the central star. These forces evolve with 
time and Fig. [S] shows them along the fiducial field line af- 
ter the stream impacts the stellar surface. At this stage, 
the pressure gradient is still effective in pushing the disk 
material out of the disk, but it acts in the opposite di- 
rection (against the free-fall of matter) in most of the ac- 
cretion stream. The pressure gradient becomes the domi- 
nant force close to the stellar surface, substantially braking 
(but not stopping) the accretion flow. As discussed below, 
at this stage of evolution, the stream has not yet reached 
a quasi-stationary condition. Later, when the stream sta- 
bilizes, the gravitational force dominates the stream evo- 
lution and the plasma continuously accelerates toward the 
star, approaching the free-fall speed tig at the star's pho- 
tosphere. Other forces acting against the free-fall of matter 
are the centrifugal force and a backward force due to the 
magnetic mi rror effect when the material approaches th e 
star (see also lRomanova et al . 2002; Z anm fc Ferreira|[2009l ). 
However, these forces are much smaller than the others and 
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Figure 5. Close view of the accretion stream during its evolution. The figure shows slices in the (R, z) plane passing through the middle 
of the stream, reporting the distributions of density (upper panels) and pressure (lower panels). The slices encompass an angular sector 
going from the disk midplane to the rotation axis. The white lines represent sampled poloidal magnetic field lines. The arrows represent 
the poloidal flow velocities. The magenta line delimits the region with plasma < 1 (on the left of each panel). 




Figure 6. Left panel: two-dimensional section, in the (x,z) plane, of the particle number density distribution, in log scale, in units of 
10 11 cm -3 (see color table in the upper left corner of the panel) at the time t = 36 hours. The slice passes through the middle of the 
accretion stream. Selected magnetic field lines are overplotted in red; the white line represents a segment of a fiducial magnetic field 
line approximately in the middle of the accretion stream. Right panel: forces governing the stream dynamics versus the distance from 
the stellar surface along the fiducial field line in the left panel: / mr is the backward force due to the magnetic mirror effect, / g the 
gravitational force, / c f the centrifugal force, /v p the pressure gradient force, and /tot is the resultant force. Forces pushing the matter 
toward (away from) the star are in blue (red). 
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Figure 7. Profiles of particle number density (upper panel) and 
relevant velocities (lower panel) versus the distance s from the 
stellar surface along the fiducial magnetic field line shown in Fig. [6] 
at time t = 36 hours: u B tr is the velocity of matter flowing along 
the accretion stream (solid black line), c a the isothermal sound 
speed (dotted red line), the free-fall velocity (dashed green 
line), and % the Alfven speed (dot-dashed blue line). 



do not play any relevant role in the stream dynamics (F ig. [6] 
see also iRomanova et al.ll2002l ;lz anm fc Ferreirallioog ). 

Figure [7] shows the profiles of particle number density 
and different velocities along the fiducial magnetic field line 
shown in Fig. [6] at time t = 36 hours. The matter flows with 
poloidal velocity u stI . The flow is accelerated by gravity and 
becomes supersonic at a distance of ~ 2 x f 11 cm from the 
disk, while the stream density decreases. The flow gradually 
approaches the free-fall velocity tiff, reaching a maximum 
velocity of u Btr ~ 200 km s _1 at s ~ fO 11 cm (u Btr ~ 0.8«ff). 
Then the flow slightly brakes while the stream density in- 
creases again approaching the stellar surface. As discussed 
above, the slowdown of the flow is due to the pressure gradi- 
ent force acting against the free-fall of matter for s ranging 
between 0.5 x fO 11 and f x 10 11 cm (see Fig. [6). This feature 
is present until the end of our 3D simulation at t ~ 48 hours. 

To further investigate the evolution of the accretion 
stream, we performed an additional simulation with the 
same parameters of the 3D simulation discussed above but 
carried out in 2.5 dimensions (2.5D), that is, in spherical co- 
ordinates (R, 6) assuming axisymmetry around the rotation 
axis of the star-disk system. Note that, in this configuration, 
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Figure 8. Profiles of velocity along the stream in the 2.5 D sim- 
ulation. As in Fig. [7] for the profiles of relevant velocities with the 
distance s from the stellar surface along a fiducial magnetic field 
line nested within the stream at time t = 36 hours (upper panel) 
and t = 68 hours (lower panel). 



the heat pulse is not localized in a relatively small portion 
of the disk (as in our 3D simulation) , but is distributed in a 
ring. Nevertheless, we found that the evolution of the flare 
and the stream described by the 2.5D simulation is quite 
similar to that of our 3D simulation. The 2.5D simulation al- 
lowed us to extend our analysis of the stream dynamics, fol- 
lowing its evolution until t = fOO hours (i.e. approximately 
four days) . Figure [8] shows the profiles of the relevant ve- 
locities derived from the 2.5D simulation at t = 36 hours 
(upper panel) and t = 68 hours (lower panel). The velocity 
profiles at t — 36 hours resemble those derived from the 3D 
simulation (lower panel in Fig. 0. In particular, the flow 
slightly brakes approaching the stellar surface due to a pres- 
sure gradient force slowing down the free-fall of matter, as 
found in the 3D simulation. In addition, the 2.5D simulation 
shows that, at this stage, the stream has not reached yet a 
quasi-stationary condition. Later, the stream stabilizes (af- 
ter t « 60 hours); the gravitational force becomes dominant 
along the stream and the matter is accelerated until it im- 
pacts the stellar surface, reaching there a maximum velocity 



of u s 



300 km s corresponding to ~ 0.9 ug. 



The above results, therefore, show that the physical 
characteristics of the accretion stream triggered by the flare 
closely recall those, largely discussed in the literature, of 
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Figure 9. Particle number density distribution on the stellar sur- 
face, in log scale, in units of 10 11 cm -3 (see color table on the 
right of the panel) at the time t = 36 hours. Selected magnetic 
field lines are overplotted in red. The white arrow marks the ro- 
tation axis of the star. 



streams driven by the accumulation of mass at the disk trun- 
cation radius under the effect of the viscosity and pushed 
out of the equa torial plane because of the g r owing pressure 
gradient there dRomanova et al.ll2002l , 120051 ; iBessolaz et al.l 
2008: IZanni fc Ferreiral20Qg| ). Our simulation shows that the 
stream is relatively cold (its temperature remains below 1 
MK). After the disk material enters the stream, its density 
slightly decreases and, close to the stellar surface, increases 
again as a result of the gas compression by the dipolar mag- 
netic field. The stream velocity u s tr gradually increases to- 
ward the star: it becomes supersonic already at a distance 
of ss 2 x 10 11 cm from the disk, and it str approaches the free- 
fall speed Uff close to the stellar surface. Figure [U] shows the 
density distribution on the stellar surface at t — 36 hours; 
the dense spot on the surface is the region of impact of the 
stream. The spot covers a small percentage of the stellar sur- 
face and the stream is inhomogeneous with its mass density 
varying across the s tream and being the larg est in the inner 
region, according to iRomanova et al.l l|2004 ). 

Finally, we derived the mass accretion rate M due to 
the stream from the side of the disk opposite to the post- 
flare loop and found M >2.5 x 10~ 10 M yr -1 . Accretion 
rates derived from optical-near UV continuum emission typ- 
ically range between 10 -12 and 10 -6 Mq yr" 1 with M vary- 
ing in the same star even by a fact or of 10 and depend- 
ing on the age and mass of the star dMuzerolle et al.l 120051 ; 
iMohantv et alj 120051 ; iNatta et akl l2006t ). In general, more 
evolved low-mass y oung stars are characterized by lower 
accretion rates (e.g. lHerczeg~fc Hillenbrand 2008 and refer- 
ences therein). We compared the rate M derived from our 3D 
simulation with the rates derived from optical-UV observa- 
tions and available in the literature. In particular, we consid- 
ered a sample of low-mass stars and brown dwarfs observed 
with LRIS on Keck I and a sample of solar-mass young ac- 
cretors observed with HST S TIS (both samples analyzed by 
iHerczeg fc Hillenbrand! [20081 . and an X-ray selected sam- 
ple of C TTSs observed with various optical telescope s (ana - 
lyzed bv lCurran et~aT1l201ll ). IHerczeg fc Hillenbrand! (^OOSt ) 
found that the excess UV and optical emission arising in the 



+ + T Tau 



+ 

+ + RU Lup 




20 



25 30 35 40 
time [ hours ] 



Figure 10. The solid line shows the mass accretion rate dur- 
ing the evolution of the funnel flow developed on the side of the 
disk opposite to the flaring loop. The crosses report the values 
of mass accretion rate derived from optical-near UV observa- 
tions for a sample of low-m ass stars and brown dwarfs (green; 
IHerczeg fc Hillenbrand! I2008T) . for a sample o f solar-mass young 
accretors (blue:~ lHerczeg fc Hill ejibranc|[200^ ) and for an X-ray- 



selected sample of CTTSs (red; ICurran et al.ll201ll) 



Balmer and Paschen continua yields mass accretion rates 
ranging between 2 x 10~ 12 and 10~ 8 Mq yr -1 in the case 
of low-mass stars and brown dwarfs, and ranging between 
2 x 1 0~ 10 and 5 x 10~ 8 M q yr -1 in the case of solar-mass 
stars. ICurran et al.l (|201 ll ) calculated M for the stars of 
their sample by measuring Ha, H/3, H 7 , Hell (4686 A), 
He I (5016 A), He I (5876 A) OI (6300 A), and He I (6678 A) 
equivalent widths and found that the mass accretion rates 
range between 2 x 10" 10 and 5 x 10" 8 Mq yr" 1 . Figure ITUI 
shows the rate M derived from our 3D simulation compared 
with the three samples of young accreting stars. The ac- 
cretion rate of our simulation is in the range of the rates 
measured in low-mass stars and lower than those of fast- 
accreting objects as BP Tau, RU Lup, or T Tau. 



4 SUMMARY AND CONCLUSIONS 

We investigated the evolution of a bright flare occurring 
close to the disk surrounding a magnetized CTTS and its 
effects on the stability of the disk, through numerical MHD 
simulations. To our knowledge, the simulation presented 
here represent the first attempt to model the 3D global evo- 
lution of the star-disk system that simultaneously include 
stellar gravity, viscosity of the disk, magnetic field, radiative 
cooling, and magnetic-field-oriented thermal conduction (in- 
cluding the effects of heat flux saturation). Our findings lead 
to several conclusions: 
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(i) During its first stage of evolution, the flare gives rise 
to a hot magnetic loop linking the inner part of the disk to 
the star; we found that the length, maximum temperature, 
and peak X-ray luminosity of the simulated flaring loop are 
similar to those derived from the analys is of luminous X-ra y 
flares observed in young low-mass stars ijFavata et al.f 2005). 

(ii) During the flare evolution, disk material evaporates in 
the outer stellar atmosphere under the effect of the thermal 
conduction. A small fraction of the evaporated disk material 
gradually fills the loop in ~ f hours. Indeed most part of 
the evaporated disk material is not efficiently confined by 
the magnetic field and channeled into the loop but is rather 
ejected away from the central star in the outer stellar corona, 
carrying away mass and angular momentum. 

(iii) In the aftermath of the flare, the disk is strongly 
perturbed: the injected heat pulse produces an overpressure 
in the disk at the loop's footpoint that travels through the 
disk. When the overpressure reaches the opposite side of the 
disk, a funnel flow starts to develop there, flowing along the 
dipolar magnetic field lines and impacting onto the stellar 
surface ~ 1 day after the injection of the heat pulse. 

(iv) We found that the mass accretion rate M of the 
stream triggered by the flare is in good agreement with 
those measured in low-mass st ars and brown dwarfs and in 
some solar-mass ac cretors fe.g. lHerczeg fe Hillenbrandll2008l : 
ICurran et al]|201ll ). The stream therefore accretes substan- 
tial mass (comparable with those observed) onto the young 
star from the side of the disk opposite to the post-flare loop. 

We conclu de therefore that the brightest flares detected 
in CTTSs (e.g. iFavata et alj 12005) can be a mechanism to 
trigger mass accretion episodes onto protostars with accre- 
tion rates in the range of those measured in low-mass stars 
and lower than those of fast-accreting objects. On the other 
hand, it is worth mentioning that large flares do not occur 
continuously in CTTSs and therefore cannot explain alone 
the time-averaged accretion rates derived in young accretors: 
an average of 1 flare per star per 650 ks (w 7 da ys) has been 
inferr ed by analyzing the COUP observ ations (Wolk et al.l 
120051 ) and clusters much older than Orion IjWolk et al.l l2004). 
To ascertain the contribution that flares may have to the 
observed time-averaged accretion rates, it is necessary to 
determine the frequency of those flares able to trigger accre- 
tion streams. In the future, an exploration of the parameter 
space of our model is therefore needed to determine to which 
extent the total energy released during the flare can be re- 
duced to still produce an accretion stream. This point is 
rather important to ascertain if a storm of small-to-medium 
flares (very frequent in CTTSs and often not even resolved 
in the lightcurves) occurring on the accretion disk of young 
stars is able to trigger accretion streams with high cadence, 
thus leading to a significant and persistent mass accretion. 
This issue deserves further investigation in future studies. 

The initial conditions adopt ed in our simulations are 
largely used in the literature (e.g. iRomanova et al.ll2002l and 
subsequent works) and are appropriate to describe the star- 
disk system in CTTSs. However, as discussed in Sect. 13.11 
these conditions may result to be idealized and simplified in 
some aspects. Many observations indicate that the stellar at- 
mospheres are permeated by magnetic fields with a degree of 
complexit y much higher than the magnetic dipole adopted 
here (e.g. iGregorv et alj lipid and references therein); the 



flare may occur in the presence of existing accretion streams 
induced by the coupling of the magnetic field to the d i sk be- 
low the corotation radius (e.g . Romanova et al.1 2002. 120031 : 
iBessolaz et al.|[2008l : IZanni fc Ferreirall2009l ) or triggered by 
other flares. In particular, the magnetic field in proximity of 
the heat release (due to magnetic reconnection) is expected 
to be more complex than that modeled here, thus confining 
more efficiently the hot plasma and enhancing the effects 
of the flare energy deposition on the system. On the other 
hand, the flaring loop and the flare-induced accretion stream 
could be strongly modified or even disappear if the flare oc- 
curs in a flux tube in itially filled with cold infalling gas. 
However, iRealel { 20091 ) showed that a flare involving dense 
infalling plasma would decay much faster than generally ob- 
served, so that the flares we address here do not occur in 
such a configuration. Nevertheless, the flaring loop and the 
flare-induced accretion stream may be strongly perturbed 
and have a more complex structure than that described here 
if the flare energy release occurs in proximity of an existing 
accretion stream. 

Following other work s in the literature (e.g. 
Romanova et al.l |2002| . 120031 . 12004 iKulkarni fc Romanova! 
20081 ) . we assumed the turbulent diffusivity to be uni- 
form within the disk, setting a as a constant. However 
numerical simulations indicate that the MRI stress would 
be substantial in the disk midplane where the plasma /3 
is high, but would drop in outer disk layers where f3 is 
relatively low. As a result, the turbulent diffusivity, scaled 
appropriately to the MRI, is expected to be smaller in the 
outer disk layers than in the disk midplane. Since the a 
value adopted in our simulation (namely a = 0.02) lies in 
the lower limit of the range of the parameter a (namel y 
0.01 < a s: 0.6; lBalbusll2003l ; IKulkarni fc Romanova! 120081 ). 
we expect that the turbulent diffusivity would be larger in 
the disk midplane than assumed in our simulation, thus 
influencing the mass accretion from the disk onto the star. 
However, the time lapse covered by our simulation (~ 48 
hours) is ~ 1/5 of the star rotation period which is much 
shorter than the time required by the tur bulent diffusivity 

(nam ely sev eral star rota t ion periods; IRomanova et all 

120021 . 120031 : IBessolaz et~ai] l200Sl : IZanni fc Ferreiral l2009h 
to be effective and trigger magnetospheric funnel streams. 
In fact, throughout the simulation, no stream develops in 
disk regions not influenced by the flare although our model 
includes the disk viscosity. 

As an additional check, we performed 2.5D simulations 
with the same parameters of the 3D simulation but exploring 
the effect of a in the range between 0.02 and 0.1, and consid- 
ering simulations either with or without the disk viscosity. 
We found that the viscosity does not significantly influence 
the results over the time laps covered by the 3D simulation. 
We conclude therefore that a bright flare occurring close 
to the circumstellar disk can trigger a substantial and long- 
lasting accretion flow with physical characteristics similar to 
those caused by the disk viscositjjf]. This new mechanism to 
drive mass accretion might be general. We suggest that the 
flaring activity common to CTTSs may turn out to be an 



2 Note that, as discussed above, IRealel {2009) has shown that 
the opposite is unlikely, that is flares cannot be triggered in an 
accreting flux tube. 
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important factor in the exchange of angular momentum and 
mass between the circumstellar disk and the central proto- 
star. This exchange is a rather hot and controversial topic 
in astrophysics and our research contributes to this issue, 
by exploring the ro le of the bright flares detected in CTTSs 
ijFavata et al.ll2005l ) not considered up to now. 

Our model predicts that the stream triggered by a flare 
impacts the stellar surface with a delay of ss 1 day since 
the peak X-ray luminosity of the flare. This time delay de- 
pends on the physical parameters of the model, namely the 
disk truncation radius, the mass density and thickness of the 
disk, and the energy released by the heat pulse. We expect 
therefore that the accretion shock generated by the impact 
of the stream on the stellar surface might be visible with a 
delay of the order of few days after the X-ray flare. Unfor- 
tunately, even if this is the case, revealing both the X-ray 
flare and the delayed accretion shock in real systems would 
be a rather difficult task because: 1) the star needs to be 
almost edge-on to allow both the flaring loop and the ac- 
cretion shock (located in different hemispheres of the star) 
to be observable; 2) both the flaring loop and the accre- 
tion shock should not be eclipsed by the star (or strongly 
absorbed by the circumstellar material) during the observa- 
tion. A more feasible task could be to detect the effects of 
the X-ray irradiation of the disk by the flare located just 
above it. For instance, a photoevaporation of disk material 
by the X-ray emission from the flare is expected to be almost 
instantaneous with the flare at variance with the case of a 
flare occurring on the star which is assumed to be delayed 
given the distance of the inner edge of the disk from the 
star. Searching for a tight correlation between X-ray flares 
and variability in IR or optical band associated with the 
photoevaporation of the disk could be an indication that 
the flares were located close to the disk. 

In addition to the results presented above, our simula- 
tion opens a number of interesting issues. We have shown 
that the flaring loop linking the disk with the star is a bright 
X-ray source, irradiating the disk from above, at variance 
with the stellar corona which irradiates only the innermost 
part of the disk. We suggest that this kind of sources may 
strongly influence the chemical and physical evolution of the 
disk, reaching parts of it that are effectively shielded from 
the stellar coronal X-ray emission, with important conse- 
quences, for instance, on the formation of planets. The tur- 
bulence and ionization state of plasma resulting from the 
irradiance of the disk by the flaring loop may also interact 
(because of the change of ionization and electrical conductiv- 
ity) in a nontrivial way with the MRI. Thus it is reasonable 
to suspect that the effects of this X-ray radiation may in- 
fluence the efficiency of MRI and, therefore, the efficiency 
of angular momentum transport within the disk. All these 
issues deserve further investigations. 
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